\(\int \frac {\cosh ^3(a+b x^2)}{x^2} \, dx\) [20]

Optimal result

Integrand size = 14, antiderivative size = 136 \[ \int \frac {\cosh ^3\left (a+b x^2\right )}{x^2} \, dx=-\frac {\cosh ^3\left (a+b x^2\right )}{x}-\frac {3}{8} \sqrt {b} e^{-a} \sqrt {\pi } \text {erf}\left (\sqrt {b} x\right )-\frac {1}{8} \sqrt {b} e^{-3 a} \sqrt {3 \pi } \text {erf}\left (\sqrt {3} \sqrt {b} x\right )+\frac {3}{8} \sqrt {b} e^a \sqrt {\pi } \text {erfi}\left (\sqrt {b} x\right )+\frac {1}{8} \sqrt {b} e^{3 a} \sqrt {3 \pi } \text {erfi}\left (\sqrt {3} \sqrt {b} x\right ) \]

[Out]

-cosh(b*x^2+a)^3/x-3/8*erf(x*b^(1/2))*b^(1/2)*Pi^(1/2)/exp(a)+3/8*exp(a)*erfi(x*b^(1/2))*b^(1/2)*Pi^(1/2)-1/8*
erf(x*3^(1/2)*b^(1/2))*b^(1/2)*3^(1/2)*Pi^(1/2)/exp(3*a)+1/8*exp(3*a)*erfi(x*3^(1/2)*b^(1/2))*b^(1/2)*3^(1/2)*
Pi^(1/2)

Rubi [A] (verified)

Time = 0.09 (sec) , antiderivative size = 136, normalized size of antiderivative = 1.00, number of steps used = 9, number of rules used = 5, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.357, Rules used = {5439, 5737, 5406, 2235, 2236} \[ \int \frac {\cosh ^3\left (a+b x^2\right )}{x^2} \, dx=-\frac {3}{8} \sqrt {\pi } e^{-a} \sqrt {b} \text {erf}\left (\sqrt {b} x\right )-\frac {1}{8} \sqrt {3 \pi } e^{-3 a} \sqrt {b} \text {erf}\left (\sqrt {3} \sqrt {b} x\right )+\frac {3}{8} \sqrt {\pi } e^a \sqrt {b} \text {erfi}\left (\sqrt {b} x\right )+\frac {1}{8} \sqrt {3 \pi } e^{3 a} \sqrt {b} \text {erfi}\left (\sqrt {3} \sqrt {b} x\right )-\frac {\cosh ^3\left (a+b x^2\right )}{x} \]

[In]

Int[Cosh[a + b*x^2]^3/x^2,x]

[Out]

-(Cosh[a + b*x^2]^3/x) - (3*Sqrt[b]*Sqrt[Pi]*Erf[Sqrt[b]*x])/(8*E^a) - (Sqrt[b]*Sqrt[3*Pi]*Erf[Sqrt[3]*Sqrt[b]
*x])/(8*E^(3*a)) + (3*Sqrt[b]*E^a*Sqrt[Pi]*Erfi[Sqrt[b]*x])/8 + (Sqrt[b]*E^(3*a)*Sqrt[3*Pi]*Erfi[Sqrt[3]*Sqrt[
b]*x])/8

Rule 2235

Int[(F_)^((a_.) + (b_.)*((c_.) + (d_.)*(x_))^2), x_Symbol] :> Simp[F^a*Sqrt[Pi]*(Erfi[(c + d*x)*Rt[b*Log[F], 2
]]/(2*d*Rt[b*Log[F], 2])), x] /; FreeQ[{F, a, b, c, d}, x] && PosQ[b]

Rule 2236

Int[(F_)^((a_.) + (b_.)*((c_.) + (d_.)*(x_))^2), x_Symbol] :> Simp[F^a*Sqrt[Pi]*(Erf[(c + d*x)*Rt[(-b)*Log[F],
2]]/(2*d*Rt[(-b)*Log[F], 2])), x] /; FreeQ[{F, a, b, c, d}, x] && NegQ[b]

Rule 5406

Int[Sinh[(c_.) + (d_.)*(x_)^(n_)], x_Symbol] :> Dist[1/2, Int[E^(c + d*x^n), x], x] - Dist[1/2, Int[E^(-c - d*
x^n), x], x] /; FreeQ[{c, d}, x] && IGtQ[n, 1]

Rule 5439

Int[Cosh[(a_.) + (b_.)*(x_)^(n_)]^(p_)*(x_)^(m_.), x_Symbol] :> Simp[-Cosh[a + b*x^n]^p/((n - 1)*x^(n - 1)), x
] + Dist[b*n*(p/(n - 1)), Int[Cosh[a + b*x^n]^(p - 1)*Sinh[a + b*x^n], x], x] /; FreeQ[{a, b}, x] && IntegersQ
[n, p] && EqQ[m + n, 0] && GtQ[p, 1] && NeQ[n, 1]

Rule 5737

Int[Cosh[w_]^(q_.)*Sinh[v_]^(p_.), x_Symbol] :> Int[ExpandTrigReduce[Sinh[v]^p*Cosh[w]^q, x], x] /; IGtQ[p, 0]
&& IGtQ[q, 0] && ((PolynomialQ[v, x] && PolynomialQ[w, x]) || (BinomialQ[{v, w}, x] && IndependentQ[Cancel[v/
w], x]))

Rubi steps \begin{align*} \text {integral}& = -\frac {\cosh ^3\left (a+b x^2\right )}{x}+(6 b) \int \cosh ^2\left (a+b x^2\right ) \sinh \left (a+b x^2\right ) \, dx \\ & = -\frac {\cosh ^3\left (a+b x^2\right )}{x}+(6 b) \int \left (\frac {1}{4} \sinh \left (a+b x^2\right )+\frac {1}{4} \sinh \left (3 a+3 b x^2\right )\right ) \, dx \\ & = -\frac {\cosh ^3\left (a+b x^2\right )}{x}+\frac {1}{2} (3 b) \int \sinh \left (a+b x^2\right ) \, dx+\frac {1}{2} (3 b) \int \sinh \left (3 a+3 b x^2\right ) \, dx \\ & = -\frac {\cosh ^3\left (a+b x^2\right )}{x}-\frac {1}{4} (3 b) \int e^{-3 a-3 b x^2} \, dx-\frac {1}{4} (3 b) \int e^{-a-b x^2} \, dx+\frac {1}{4} (3 b) \int e^{a+b x^2} \, dx+\frac {1}{4} (3 b) \int e^{3 a+3 b x^2} \, dx \\ & = -\frac {\cosh ^3\left (a+b x^2\right )}{x}-\frac {3}{8} \sqrt {b} e^{-a} \sqrt {\pi } \text {erf}\left (\sqrt {b} x\right )-\frac {1}{8} \sqrt {b} e^{-3 a} \sqrt {3 \pi } \text {erf}\left (\sqrt {3} \sqrt {b} x\right )+\frac {3}{8} \sqrt {b} e^a \sqrt {\pi } \text {erfi}\left (\sqrt {b} x\right )+\frac {1}{8} \sqrt {b} e^{3 a} \sqrt {3 \pi } \text {erfi}\left (\sqrt {3} \sqrt {b} x\right ) \\ \end{align*}

Mathematica [A] (verified)

Time = 0.25 (sec) , antiderivative size = 204, normalized size of antiderivative = 1.50 \[ \int \frac {\cosh ^3\left (a+b x^2\right )}{x^2} \, dx=\frac {-6 \cosh \left (a+b x^2\right )-2 \cosh \left (3 \left (a+b x^2\right )\right )+3 \sqrt {b} \sqrt {\pi } x \cosh (a) \text {erfi}\left (\sqrt {b} x\right )+\sqrt {b} \sqrt {3 \pi } x \cosh (3 a) \text {erfi}\left (\sqrt {3} \sqrt {b} x\right )+3 \sqrt {b} \sqrt {\pi } x \text {erfi}\left (\sqrt {b} x\right ) \sinh (a)+3 \sqrt {b} \sqrt {\pi } x \text {erf}\left (\sqrt {b} x\right ) (-\cosh (a)+\sinh (a))+\sqrt {b} \sqrt {3 \pi } x \text {erfi}\left (\sqrt {3} \sqrt {b} x\right ) \sinh (3 a)+\sqrt {b} \sqrt {3 \pi } x \text {erf}\left (\sqrt {3} \sqrt {b} x\right ) (-\cosh (3 a)+\sinh (3 a))}{8 x} \]

[In]

Integrate[Cosh[a + b*x^2]^3/x^2,x]

[Out]

(-6*Cosh[a + b*x^2] - 2*Cosh[3*(a + b*x^2)] + 3*Sqrt[b]*Sqrt[Pi]*x*Cosh[a]*Erfi[Sqrt[b]*x] + Sqrt[b]*Sqrt[3*Pi
]*x*Cosh[3*a]*Erfi[Sqrt[3]*Sqrt[b]*x] + 3*Sqrt[b]*Sqrt[Pi]*x*Erfi[Sqrt[b]*x]*Sinh[a] + 3*Sqrt[b]*Sqrt[Pi]*x*Er
f[Sqrt[b]*x]*(-Cosh[a] + Sinh[a]) + Sqrt[b]*Sqrt[3*Pi]*x*Erfi[Sqrt[3]*Sqrt[b]*x]*Sinh[3*a] + Sqrt[b]*Sqrt[3*Pi
]*x*Erf[Sqrt[3]*Sqrt[b]*x]*(-Cosh[3*a] + Sinh[3*a]))/(8*x)

Maple [A] (verified)

Time = 0.18 (sec) , antiderivative size = 149, normalized size of antiderivative = 1.10

method	result	size
risch	\(-\frac {{\mathrm e}^{-3 a} {\mathrm e}^{-3 b \,x^{2}}}{8 x}-\frac {{\mathrm e}^{-3 a} \sqrt {b}\, \sqrt {\pi }\, \sqrt {3}\, \operatorname {erf}\left (x \sqrt {3}\, \sqrt {b}\right )}{8}-\frac {3 \,{\mathrm e}^{-a} {\mathrm e}^{-b \,x^{2}}}{8 x}-\frac {3 \,\operatorname {erf}\left (x \sqrt {b}\right ) \sqrt {b}\, \sqrt {\pi }\, {\mathrm e}^{-a}}{8}-\frac {{\mathrm e}^{3 a} {\mathrm e}^{3 b \,x^{2}}}{8 x}+\frac {3 \,{\mathrm e}^{3 a} b \sqrt {\pi }\, \operatorname {erf}\left (\sqrt {-3 b}\, x \right )}{8 \sqrt {-3 b}}-\frac {3 \,{\mathrm e}^{a} {\mathrm e}^{b \,x^{2}}}{8 x}+\frac {3 \,{\mathrm e}^{a} b \sqrt {\pi }\, \operatorname {erf}\left (\sqrt {-b}\, x \right )}{8 \sqrt {-b}}\)	\(149\)

[In]

int(cosh(b*x^2+a)^3/x^2,x,method=_RETURNVERBOSE)

[Out]

-1/8/exp(a)^3/x*exp(-3*b*x^2)-1/8/exp(a)^3*b^(1/2)*Pi^(1/2)*3^(1/2)*erf(x*3^(1/2)*b^(1/2))-3/8/exp(a)/x*exp(-b
*x^2)-3/8*erf(x*b^(1/2))*b^(1/2)*Pi^(1/2)/exp(a)-1/8*exp(a)^3/x*exp(3*b*x^2)+3/8*exp(a)^3*b*Pi^(1/2)/(-3*b)^(1
/2)*erf((-3*b)^(1/2)*x)-3/8*exp(a)*exp(b*x^2)/x+3/8*exp(a)*b*Pi^(1/2)/(-b)^(1/2)*erf((-b)^(1/2)*x)

Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 891 vs. \(2 (98) = 196\).

Time = 0.26 (sec) , antiderivative size = 891, normalized size of antiderivative = 6.55 \[ \int \frac {\cosh ^3\left (a+b x^2\right )}{x^2} \, dx=\text {Too large to display} \]

[In]

integrate(cosh(b*x^2+a)^3/x^2,x, algorithm="fricas")

[Out]

-1/8*(cosh(b*x^2 + a)^6 + 6*cosh(b*x^2 + a)*sinh(b*x^2 + a)^5 + sinh(b*x^2 + a)^6 + 3*(5*cosh(b*x^2 + a)^2 + 1
)*sinh(b*x^2 + a)^4 + 3*cosh(b*x^2 + a)^4 + 4*(5*cosh(b*x^2 + a)^3 + 3*cosh(b*x^2 + a))*sinh(b*x^2 + a)^3 + sq
rt(3)*sqrt(pi)*(x*cosh(b*x^2 + a)^3*cosh(3*a) + x*cosh(b*x^2 + a)^3*sinh(3*a) + (x*cosh(3*a) + x*sinh(3*a))*si
nh(b*x^2 + a)^3 + 3*(x*cosh(b*x^2 + a)*cosh(3*a) + x*cosh(b*x^2 + a)*sinh(3*a))*sinh(b*x^2 + a)^2 + 3*(x*cosh(
b*x^2 + a)^2*cosh(3*a) + x*cosh(b*x^2 + a)^2*sinh(3*a))*sinh(b*x^2 + a))*sqrt(-b)*erf(sqrt(3)*sqrt(-b)*x) + sq
rt(3)*sqrt(pi)*(x*cosh(b*x^2 + a)^3*cosh(3*a) - x*cosh(b*x^2 + a)^3*sinh(3*a) + (x*cosh(3*a) - x*sinh(3*a))*si
nh(b*x^2 + a)^3 + 3*(x*cosh(b*x^2 + a)*cosh(3*a) - x*cosh(b*x^2 + a)*sinh(3*a))*sinh(b*x^2 + a)^2 + 3*(x*cosh(
b*x^2 + a)^2*cosh(3*a) - x*cosh(b*x^2 + a)^2*sinh(3*a))*sinh(b*x^2 + a))*sqrt(b)*erf(sqrt(3)*sqrt(b)*x) + 3*sq
rt(pi)*(x*cosh(b*x^2 + a)^3*cosh(a) + x*cosh(b*x^2 + a)^3*sinh(a) + (x*cosh(a) + x*sinh(a))*sinh(b*x^2 + a)^3
+ 3*(x*cosh(b*x^2 + a)*cosh(a) + x*cosh(b*x^2 + a)*sinh(a))*sinh(b*x^2 + a)^2 + 3*(x*cosh(b*x^2 + a)^2*cosh(a)
+ x*cosh(b*x^2 + a)^2*sinh(a))*sinh(b*x^2 + a))*sqrt(-b)*erf(sqrt(-b)*x) + 3*sqrt(pi)*(x*cosh(b*x^2 + a)^3*co
sh(a) - x*cosh(b*x^2 + a)^3*sinh(a) + (x*cosh(a) - x*sinh(a))*sinh(b*x^2 + a)^3 + 3*(x*cosh(b*x^2 + a)*cosh(a)
- x*cosh(b*x^2 + a)*sinh(a))*sinh(b*x^2 + a)^2 + 3*(x*cosh(b*x^2 + a)^2*cosh(a) - x*cosh(b*x^2 + a)^2*sinh(a)
)*sinh(b*x^2 + a))*sqrt(b)*erf(sqrt(b)*x) + 3*(5*cosh(b*x^2 + a)^4 + 6*cosh(b*x^2 + a)^2 + 1)*sinh(b*x^2 + a)^
2 + 3*cosh(b*x^2 + a)^2 + 6*(cosh(b*x^2 + a)^5 + 2*cosh(b*x^2 + a)^3 + cosh(b*x^2 + a))*sinh(b*x^2 + a) + 1)/(
x*cosh(b*x^2 + a)^3 + 3*x*cosh(b*x^2 + a)^2*sinh(b*x^2 + a) + 3*x*cosh(b*x^2 + a)*sinh(b*x^2 + a)^2 + x*sinh(b
*x^2 + a)^3)

Sympy [F]

\[ \int \frac {\cosh ^3\left (a+b x^2\right )}{x^2} \, dx=\int \frac {\cosh ^{3}{\left (a + b x^{2} \right )}}{x^{2}}\, dx \]

[In]

integrate(cosh(b*x**2+a)**3/x**2,x)

[Out]

Integral(cosh(a + b*x**2)**3/x**2, x)

Maxima [A] (verification not implemented)

none

Time = 0.26 (sec) , antiderivative size = 102, normalized size of antiderivative = 0.75 \[ \int \frac {\cosh ^3\left (a+b x^2\right )}{x^2} \, dx=-\frac {\sqrt {3} \sqrt {b x^{2}} e^{\left (-3 \, a\right )} \Gamma \left (-\frac {1}{2}, 3 \, b x^{2}\right )}{16 \, x} - \frac {\sqrt {3} \sqrt {-b x^{2}} e^{\left (3 \, a\right )} \Gamma \left (-\frac {1}{2}, -3 \, b x^{2}\right )}{16 \, x} - \frac {3 \, \sqrt {b x^{2}} e^{\left (-a\right )} \Gamma \left (-\frac {1}{2}, b x^{2}\right )}{16 \, x} - \frac {3 \, \sqrt {-b x^{2}} e^{a} \Gamma \left (-\frac {1}{2}, -b x^{2}\right )}{16 \, x} \]

[In]

integrate(cosh(b*x^2+a)^3/x^2,x, algorithm="maxima")

[Out]

-1/16*sqrt(3)*sqrt(b*x^2)*e^(-3*a)*gamma(-1/2, 3*b*x^2)/x - 1/16*sqrt(3)*sqrt(-b*x^2)*e^(3*a)*gamma(-1/2, -3*b
*x^2)/x - 3/16*sqrt(b*x^2)*e^(-a)*gamma(-1/2, b*x^2)/x - 3/16*sqrt(-b*x^2)*e^a*gamma(-1/2, -b*x^2)/x

Giac [F]

\[ \int \frac {\cosh ^3\left (a+b x^2\right )}{x^2} \, dx=\int { \frac {\cosh \left (b x^{2} + a\right )^{3}}{x^{2}} \,d x } \]

[In]

integrate(cosh(b*x^2+a)^3/x^2,x, algorithm="giac")

[Out]

integrate(cosh(b*x^2 + a)^3/x^2, x)

Mupad [F(-1)]

Timed out. \[ \int \frac {\cosh ^3\left (a+b x^2\right )}{x^2} \, dx=\int \frac {{\mathrm {cosh}\left (b\,x^2+a\right )}^3}{x^2} \,d x \]

[In]

int(cosh(a + b*x^2)^3/x^2,x)

[Out]

int(cosh(a + b*x^2)^3/x^2, x)